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Abstract — the aim of this paper is on the problem of crack 
detection for a cracked beam using theoretical and software 
analysis. Crack is the discontinuation in a body. Location of 
crack through the method of vibration is great advantages over 
traditional methods. Presence of crack in a structural member 
induces local flexibility which affects vibration response of the 
structure. This property may be used to detect location and 
depth of the crack in structural member. The presence of crack 
leads to changes in some of the lower natural frequencies and 
mode shapes. Method based on measurement of natural 
frequencies is presented for detection of the location and size of 
a crack in a cantilever beam. Numerical calculations has been 
done by solving the Euler equation for un-crack beam and 
cracked beam to obtain first three natural frequencies of 
different modes of vibration considering various crack 
positions for the beam. CATIA software is used for designing 
cantilever beam. ANSYS software is used for analysis of crack 
and un-crack cantilever beam. Total 7 models has been tested 
according to theoretical and software method. 

Index terms — analysis, Finite element analysis, structure. 
I. I. INTRODUCTION 


D. Specifications: 


Cantilever beam height- 10mm 

Crack depth in cantilever beam - 1mm, 2mm, and 3mm 
Crack location- in cantilever beam 150mm, 300mm 
E = modulus of elasticity (N/m 2 ) 
m = mass of the beam per unit length (kg/m) 

1= moment of inertia (m 4 ). 

coi = natural frequency of the i th mode (rad/s) 

Cantilever beam length- 600mm 
Cantilever beam width- 10mm 
Ll=Distance between the left node and crack, 
I 0 =Moment of inertia of the beam cross section, 

I c = Moment of inertia of the beam with crack. 
q=displacement vector of the beam 



is the Stiffness matrix of the cracked element 


| K e| 

L - 1 is the Element stiffness matrix, 
L= length of the beam, 


[KcJ 


is the Reduction in stiffness matrix due to the crack. 


A. Vibration: 

Motion which repeats itself after a certain interval of time is 
called vibration. A vibration can caused due to external 
unbalanced force also. A vibratory system, in general, 
includes elastic member for storing potential 

B. Lateral vibration: 

The vibration in which particles of the system vibrates in the 
direction perpendicular to the axis of system is known as 
lateral vibration. Since the cantilever beam has infinite 
number of masses, we need infinite number of coordinates to 
specify the deflected configuration. Thus, cantilever beam is 
infinite degree of freedom system and it is necessary to study 
lateral vibration of cantilever beam. 

C. Need of Crack Detection: 

Crack detection is one of the important aspects in structural 
engineering for safety reasons and because of economic 
benefits that can result. 
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II. LITERATURE REVIEW 

Vibration analysis of a cracked beam was performed by S 
Orhan et al. in order to identify the crack in a cantilever 
beam. Single- and two-edge cracks were evaluated. The study 
results suggest that free vibration analysis provides suitable 
information for the detection of single and two cracks, 
whereas forced vibration can detect only the single crack 
condition. The Euler- Bernoulli beam model was assumed. 
The crack is assumed to be an open crack and the damping 
has not been considered in this study. 

E Cam, S Orhan and M Luy proposed method of analysis of 
cracked beam structure using impact echo method. They 
analyzed vibrations as a result of impact shocks. ANSYS 
software was used for experimental results and simulations. 
R. Tiwari and M. Karthikey developed identification 
procedure for the detection, localization, and size of a crack 
in a beam based on forced response measurements. They used 
circular beam supported by rolling bearings at both ends for 
this experiment. 

C. P. Filipich, M B Rosales and F S Buezas et al. determined 
sensibility analysis of the inverse problem of the crack 
parameters (location and depth). M. Jabalpur, N. Khaji, M. 
Shafiei developed analytical approach for crack 
identification procedure in uniform beams with an open edge 
crack, based on bending vibration measurements. S. 
Loutridis, E.Douk developed new method for crack detection 
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in beams based on instantaneous frequency and empirical 
mode decomposition. It follows that the harmonic distortion 
increases with crack depth following definite trends and can 
be also used as an effective indicator for crack size. The 
research work by W Zhang, Z Wang and H Ma et al. 
illustrates the crack identification method combining wavelet 
analysis with transform matrix. 

III. THEROTICAL FORMULATIONS 

Dynamic response of the structure affected by the following 
aspects of the crack 

• Position of crack 

• Depth of crack 

• Number of cracks 

A. Equation of free vibration of beam: 

The beam with a transverse edge crack is clamped at left 
end, free at right end and has different cross section The 
crack in this case is assumed to be an open crack and the 
damping is not being considered in this theory. 


The element stiffness matrix of the un-cracked beam is given 
as: 

[K e ]= i[B(x)] T EI[B(x)]dx 


[B(x)]=(H,(x)H 2 (x)H 3 (x)H 4 (x)} 


(4) 


Where the Hermitian shape functions are defined as: 

3x 2 2x 3 

Hj (x) = 1 - — h 

\Z 


H 2 (x) = 




3x 2 2x 3 

H 3 (x) = -^— J- 


H 4 (x) = - 
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Figure A: Cantilever Beam without Crack. 
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Figure B: Cantilever Beam with Crack. 


The free bending vibration of an Euler-Bernoulli beam of a 
constant rectangular cross section is given by the following 
equation: 


EI-Up-mco 2 y = 0 
dx 4 


( 1 ) 


i 4 _ mto f 

By defining EI equation is rearranged as a 

fourth-order differential equation as follows: 


Assuming the beam rigidity EI is constant and is given by EI 0 
within the element, and then the element stiffness is: 



' 12 

61 

-12 

61 

El 0 

61 

41 2 

-61 

21 2 

l 3 

-12 

-61 

12 

-61 


61 

21 2 

-61 

41 2 


[K*l = [K e ]-|K c | 


Here, 

The matrix [Kc] is: 
Ki 


|K C 


Ml 
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12 


-K 


Where: 
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The general solution to the equation is: 

y=AcosU+BsinU+CcoshLx+Dsinhl:x 

1111 (3) 

Where A, B, C, D are constants and 6 V is a frequency 
parameter. 
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12E(I 0 -I C ) 


K 24= t2 


4 + 21 , 


( 

9Li 9L: 
2-4 +_L 

L L 2 


21 


K44 = 


12E(I 0 -I C ) 


31c 

(3Li V 

4+21 

-U 

[l 2 

1 Jj 


Here, lc=1.5W, It is supposed that the crack does not affect 
the mass distribution of the beam. Therefore, the consistent 
mass matrix of the beam element can be formulated directly 
as: 

i 

[m c ]= f pA[H(x)] T [H(x)]dx 

0 (7) 
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The natural frequency then can be calculated from the 
relation: 


[-G) 2 [M]+[K]]{q} = 0 


(9) 


From Euler’s Beam theory theoretical formula for natural 
frequency given as follows: 


Where C is Constant, 

C = 0.56 for First mode, 3.51 for Second, 9.82 for Third 
mode for Cantilever Beam. 


IV. FINITE ELEMENT ANALYSIS 

Finite element analysis has been carried out by ANSYS 
software. 

In general, a finite-element solution may be broken into the 
following three stages. 


Following steps show the guidelines for carrying out Modal 
analysis. 

Define Materials 

1. Set preferences. (Structural) 

2. Define constant material properties. 

Model the Geometry 

3. Follow bottom up modeling and create/import the 
geometry 

Generate Mesh 

4. Define element type. 

5. Mesh the area. 

Apply Boundary Conditions 

6. Apply constraints to the model. 

Obtain Solution 

7. Specify analysis types and options. 

8. Solve. 

Total 7 models has been tested .The result of ANSYS 
analysis for the beams have first natural frequencies is shown 
below: 



Figure A: Result of beam model 1 



Figure B: Result of beam model 2 


A. Pre-processing: 

The major steps in pre-processing are 

(i) Define key points/lines/areas/volumes, 

(ii) Define element type and material properties and (iii) 
mesh lines/areas/ volumes as required. 

B. Solution : 

Assigning loads, constraints and solving. Here, it is 
necessary to specify the loads, constraints and finally solve 
the resulting set of equation. 

C. Post processing : 

In this stage we can see (i) lists of nodal displacements, (ii) 
element forces and moments, (iii) deflection plots and (iv) 
frequencies and temperature maps. 



Figure C: Result of beam model 3 
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V. RESULTS 

The table 1 shows theoretical Natural Frequency, Crack 

Depth and Location 

Where FNF = first natural frequency. 

SNF = second natural frequency. 

TNF= third natural frequency. 


Table 1 


Sr. 

No. 

Location 
of crack 
(mm) 

Depth 

of 

crack 

(mm) 

FNF 

(Hz) 

SNF 

(Hz) 

TNF 

(Hz) 

1 

0 

0 

22.661 

142.441 

397.376 

2 

150 

1 

22.652 

142.433 

397.375 

3 

150 

2 

22.566 

141.839 

395.712 

4 

150 

3 

22.350 

40.506 

391.981 

5 

300 

1 

22.651 

142.433 

397.374 

6 

300 

2 

22.566 

141.837 

395.715 

7 

300 

3 

22.353 

40.506 

391.981 


The table 2 shows results of beams tested by finite element 
modeling. 


Table 2 


Sr. 

No. 

Location 
of crack 
(mm) 

Depth 
of crack 
(mm) 

FNF 

(Hz) 

SNF 

(Hz) 

TNF 

(Hz) 

1 

0 

0 

22.661 

141.84 

396.32 

2 

150 

1 

22.648 

141.84 

395.82 

3 

150 

2 

22.613 

141.82 

394.41 

4 

150 

3 

22.561 

141.77 

395.16 

5 

300 

1 

22.641 

141.75 

396.34 

6 

300 

2 

22.623 

140.78 

396.33 

7 

300 

3 

22.572 

141.17 

396.33 
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Graph A: First Natural Frequency Vs Crack Depth 
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We plot the Graph A. 

Where X axis = Crack Depth (mm) 

Y axis= First Natural Frequency (Hz) 

VI. CONCLUSION 

From Table 1 , Table 2 and Graph A: 

• Natural frequency of the beam decreases as the crack depth 

increases and the crack location is constant. 

• Changes in natural frequency observed at the area of crack 

location 

• The results show that the values of natural frequencies by 

theory and ANSYS are close to the agreement. 

• Natural frequencies vary as the crack location from the 

cantilever end varied 
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